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Abstract First, we show how the quantum circuits for generating and measuring multi-party 
entanglement of qubits can be translated to continuous quantum variables. We 
derive sufficient inseparability criteria for iV-party continuous-variable states and 
discuss their applicability. Then, we consider a family of multipartite entangled 
states (multi-party multi-mode states with one mode per party) described by 
continuous quantum variables and analyze their properties. These states can be 
efficiently generated using squeezed light and linear optics. 
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1. INTRODUCTION 

What is the main motivation to deal with continuous variables for quantum 
communication purposes? Quantum communication schemes rely on state 
preparation, local unitary transformations, measurements, and classical com- 
munication. In addition, sometimes shared entanglement is part of the protocol. 
Within the framework of quantum optics, these ingredients can be efficiently 
implemented when they are applied to the continuous quadrature amplitudes 
of electromagnetic modes. For example, the tools for measuring a quadrature 
with near-unit efficiency or for displacing an optical mode in phase space are 
provided by homodyne detection and feed-forward techniques, respectively. 
Continuous-variable entanglement can be efficiently produced using squeezed 
light and linear optics. In this chapter, we consider a rather general manifes- 
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tation of continuous-variable entanglement, namely that between an arbitrary 
number of parties (modes). We will see that even those N -party entangled 
states where none of the N parties can be separated from the others in the 
total state vector are comparatively "cheap" in terms of the resources needed: 
their generation only requires one single-mode squeezed state and N — 1 beam 
splitters. 

2. MULTIPARTITE ENTANGLEMENT 

The main subject of this section is multi-party entanglement of infinite- 
dimensional states described by continuous variables. After a few general 
remarks on entanglement between two and more parties in arbitrary dimen- 
sions, we will show how the quantum circuits for creating and measuring qubit 
entanglement may be translated to continuous variables. Then we derive in- 
equalities that may serve as sufficient multi-party inseparability criteria for 
continuous- variable states. These are applicable both for a theoretical test and 
for an indirect experimental verification of multi-party entanglement. Finally, 
we focus on a family of genuinely multi-party entangled continuous-variable 
states whose members are fully inseparable with respect to all their parties. 

2.1 TWO PARTIES VERSUS MANY PARTIES 

Bipartite entanglement, the entanglement of a pair of systems shared by 
two parties, is easy to handle for pure states. For any pure two-party state, 
orthonormal bases of each subsystem exist, {|u n )} an d {\v n )},so that the total 
state vector can be written in the "Schmidt decomposition" [1] as 

W) = ^2c n \u n }\v n ) , (1.1) 

n 

where the summation goes over the smaller of the dimensionalities of the two 
subsystems. The Schmidt coefficients c n are real and non-negative, and satisfy 
Y2 n c n = 1. The Schmidt decomposition may be obtained by transforming the 
expansion of an arbitrary pure bipartite state as 

= ^2a m k\m)\k) = ^2 UmnCnnVkn\m)\k) = ^2 c n \u n )\v n ) , (1.2) 

mk nmk n 

with c nn = c n . In the first step, the matrix a with complex elements a m k 
is diagonalised, a = ucv T , where u and v are unitary matrices and c is a 
diagonal matrix with non-negative elements. In the second step, we defined 

\ u n) = Em 

) and \v n ) = J2k v kn\k) which form orthonormal sets due 
to the unitarity of u and v and the orthonormality of \m) and | k) . A pure state of 
two ci-level systems ("qudits") is now maximally entangled when the Schmidt 
coefficients of its total state vector are all equal. Since the eigenvalues of the 
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reduced density operator upon tracing out one half of a bipartite state are the 
Schmidt coefficients squared, 

pi = Tr 2 pi2 = Tr 2 |^)i2(^| = ^ c l\ u n)i(u n \ , (1.3) 

n 

tracing out either qudit of a maximally entangled state leaves the other half 
in the maximally mixed state 11 jd. A pure two-party state is factorizable (not 
entangled) if and only if the number of nonzero Schmidt coefficients is one. Any 
Schmidt number greater than one indicates entanglement. Thus, the "majority" 
of pure state vectors in the Hilbert space of two parties are nonmaximally 
entangled. Furthermore, any pure two-party state is entangled if and only if 
for suitably chosen observables, it yields a violation of inequalities imposed by 
local realistic theories [2]. A unique measure of bipartite entanglement for pure 
states is given by the partial von Neumann entropy, the von Neumann entropy 
(— Trplog/5) of the remaining system after tracing out either subsystem [3]: 
-Ev.N. = — Trpi log d pi = — J2n c n l°gd c n» ranging between zero and one (in 
units of "edits"). 

Mixed states are more subtle, even for only two parties. As for the quan- 
tification of bipartite mixed-state entanglement, there are various measures 
available such as the entanglement of formation and distillation [4]. Only for 
pure states, these measures coincide and equal the partial von Neumann en- 
tropy. The definition of pure-state entanglement via the non-factorizability of 
the total state vector is generalized to mixed states through non-separability 
(or inseparability) of the total density operator. A general quantum state of a 
two-party system is separable if its total density operator is a mixture (a convex 
sum) of product states [5], 

Pl2 = ^PiPil® Pi2 ■ (1-4) 

i 

Otherwise, it is inseparable 1 . In general, it is a non-trivial question whether 
a given density operator is separable or inseparable. Nonetheless, a very 
convenient method to test for inseparability is Peres' partial transpose criterion 
[6]. For a separable state as in Eq. (1.4), transposition of either density matrix 



'Separable states also exhibit correlations, but those are purely classical. For instance, compare the 
separable state p = ^ (|0> <0| (g> |0)(0| + |1}(1| (g) |1)(1|) to the pure maximally entangled "Bell state" 
1*+} = ^tjOO} ® |0) + |1) ® |1>) = ^75(1+} ® |+) + |-> ® |-» with the conjugate basis states 
|±) = -^=(|0) ± |1)). The separable state p is classically correlated only with respect to the predetermined 

basis {|0), |1)}. However, the Bell state |$ + ) is a priori quantum correlated in both bases {|0}, |1}} and 
{|+), |— )}, and may become a posteriori classically correlated depending on the particular basis choice in a 
local measurement. Similarly, we will see later that the inseparability criteria for continuous variables need 
to be expressed in terms of the positions and their conjugate momenta. 
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yields again a legitimate non-negative density operator with unit trace, 

Pl2 = ^Piipilf ® Pi2 , (1-5) 

i 

since {pn) T = (pn)* corresponds to a legitimate density matrix. This is a nec- 
essary condition for a separable state, and hence a single negative eigenvalue 
of the partially transposed density matrix is a sufficient condition for insepa- 
rability. In the (2 x 2)- and (2 x 3)-dimensional cases (and, for example, for 
two-mode Gaussian states, see below), this condition is both necessary and suf- 
ficient. For any other dimension, negative partial transpose is only sufficient for 
inseparability [7] 2 . Other sufficient inseparability criteria include violations 
of inequalities imposed by local realistic theories (though mixed inseparable 
states do not necessarily lead to such violations), an entropic inequality [namely 
-Ev.n.(Pi) > ^v.N.(pi2)> again with p\ = TV2P12] [11], and a condition based 
on the theory of majorization [12]. Concluding the discussion of two-party 
entanglement, we emphasize that both the pure-state Schmidt decomposition 
and the partial transpose criterion for mixed states are also applicable to infinite 
dimensions. An example for the infinite-dimensional Schmidt decomposition 
is the two-mode squeezed vacuum state in the Fock (photon number) basis 
[13]. The unphysical operation (a positive, but not completely positive map) 
that corresponds to the transposition is time reversal [14]: in terms of contin- 
uous variables, any separable two-party state remains a legitimate state after 
the transformation (x\,pi,X2,P2) — ► (xi, —Pi,X2,P2), where (xi,pi) are the 
phase-space variables (positions and momenta) for example in the Wigner 
representation. However, arbitrary inseparable states may be turned into un- 
physical states, and furthermore, inseparable two-party two-mode Gaussian 
states always become unphysical via this transformation [14]. 

Multipartite entanglement, the entanglement shared by more than two 
parties, is a more complex issue. For pure multi-party states, a Schmidt 
decomposition does not exist in general. The total state vector then cannot be 
written as a single sum over orthonormal basis states. There is, however, one 
very important representative of multipartite entanglement which does have the 
form of a multi-party Schmidt decomposition, namely the Greenberger-Horne- 



2 Inseparable states with positive partial transpose cannot be distilled to a maximally entangled state via local 
operations and classical communication. They are so-called "bound entangled" [8]. The converse, however, 
does not hold. An explicit example of a bound entangled state with negative partial transpose was given 
in Ref . [9] . In other words, not all entangled states that reveal their inseparability through negative partial 
transpose are distillable or "free entangled". On the other hand, any state pi2 that violates the so-called 
reduction criterion, pi (g> fl. — pi2 > or ll Cg) p2 — pi2 > 0, is both inseparable and distillable [10]. This 
reduction criterion is in general weaker than the partial transpose criterion and the two criteria are equivalent 
in the (2 X 2)- and (2 X 3)-dimensional cases. 
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Zeilinger (GHZ) state [15] 

|GHZ) = -L(|000> + |111» , (1.6) 

here given as a three-qubit state. Although there is no rigorous definition of 
maximally entangled multi-party states due to the lack of a general Schmidt 
decomposition, the form of the GHZ state with all "Schmidt coefficients" equal 
suggests that it exhibits maximum multipartite entanglement. In fact, there 
are various reasons for assigning the attribute "maximally entangled" to the 
iV-party GHZ states [(|000 • • • 000) + | 111 • • • lll))/>/2]. For example, they 
yield the maximum violations of multi-party inequalities imposed by local 
realistic theories [16]. Further, their entanglement heavily relies on all parties, 
and if examined pairwise they do not contain simple bipartite entanglement 
(see below). 

For the case of three qubits, any pure and fully entangled state can be 
transformed to either the GHZ state or the so-called W state [17], 

|W> = ^=(1100) + |010) + |001)) , (1.7) 

via stochastic local operations and classical communication ("SLOCC", where 
stochastic means that the state is transformed with non-zero probability). Thus, 
with respect to SLOCC, there are two inequivalent classes of genuine tripartite 
entanglement, represented by the GHZ and the W state. Genuinely or fully 
tripartite entangled here means that the entanglement of the three-qubit state is 
not just present between two parties while the remaining party can be separated 
by a tensor product. Though genuinely tripartite, the entanglement of the W 
state is also "readily bipartite". This means that the remaining two-party state 
after tracing out one party, 

Tn|W)(W| = 1 (|00)(00| + |10)(10| + |01)<01| + |01)(10| + |10)<01|) , 

(1.8) 

is inseparable which can be verified by taking the partial transpose [the eigen- 
values are 1/3, 1/3, (1 ± \/5)/6]. This is in contrast to the GHZ state where 
tracing out one party yields the separable two-qubit state 

Tn|GHZ)(GHZ| = i(|00)(00| +|11)(11|) (1.9) 

= ^(|0><0|®|0><0| + |1><1|®|1><1|) . 

Note that this is not the maximally mixed state of two qubits, l® 2 /4. The 
maximally mixed state of one qubit, however, is obtained after tracing out two 
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parties of the GHZ state. Maximum bipartite entanglement is available from 
the GHZ state through a local measurement of one party in the conjugate basis 
{(|0) ± |l))/\/2} (plus classical communication about the result), 



i (I0h ±IDi) (i<0|±i<H)|GHZ) l :(|o)i±|1)i)0| ^ 



(|0) 1 ±|1) 1 )( 1 (0|± 1 (1|)|GHZ)|| V2 K 



(1.10) 



Here, |$±) are two of the four Bell states, |$±) = (|00) ± |ll))/>/2, 1^) = 
(|01) ± \10))/V2. 

What can be said about arbitrary mixed entangled states of more than two 
parties? There is of course an immense variety of inequivalent classes of multi- 
party mixed states [e.g., five classes of three-qubit states of which the extreme 
cases are the fully separable (p = J2i Pi Pn® Pi2® Pa) and the fully (genuinely) 
inseparable states [18]]. In general, multi-party inseparability criteria cannot be 
formulated in such a compact form as the two-party partial transpose criterion. 
Similarly, the quantification of multipartite entanglement, even for pure states, 
is still subject of current research. Existing multi-party entanglement measures 
do not appear to be unique as is the partial von Neumann entropy for pure 
two-party states. Furthermore, violations of multi-party inequalities imposed 
by local realism do not necessarily imply genuine multi-party inseparability. In 
the case of continuous variables, we may now focus on the following questions: 
How can we generate, measure, and (theoretically and experimentally) verify 
genuine multipartite entangled states? How do the continuous-variable states 
compare to their qubit counterparts with respect to various properties? 

2.2 CREATING MULTIPARTITE ENTANGLEMENT 

A compact way to describe how entanglement may be created is in terms 
of a quantum circuit. Quantum circuits consist of a sequence of unitary trans- 
formations (quantum gates), sometimes supplemented by measurements. A 
quantum circuit is independent of a particular physical realization. 

Let us consider the generation of entanglement between arbitrarily many 
qubits. The quantum circuit shall turn N independent qubits into an iV-partite 
entangled state. Initially, the N qubits shall be in the eigenstate |0). All we 
need is a circuit with the following two elementary gates: the Hadamard gate, 
acting on a single qubit as 

|o>— >^(|0> + |1», |i> — ^(|o>-|i». (lid 

and the controlled-NOT (C-NOT) gate, a two-qubit operation acting as 
|00> — |00>, |01> — |01>, |10> — |11>, |H) — |10>. (1.12) 
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Figure 1.1 Quantum circuit for generating the iV-qubit GHZ state. The gates (unitary trans- 
formations) are a Hadamard gate ("H") and pairwise acting C-NOT gates. 



The first qubit (control qubit) remains unchanged under the C-NOT. The second 
qubit (target qubit) is flipped if the control qubit is set to 1, and is left unchanged 
otherwise. Equivalently, we can describe the action of the C-NOT gate by 
\yi,U2) — > 1 2/1 , 2/1 © 2/2) with 2/1,2/2 = 0,1 and the addition modulo two ©. 
The iV-partite entangled output state of the circuit (see Fig. 1.1) is the iV-qubit 
GHZ state. 

Let us translate the qubit quantum circuit to continuous variables [19]. 
The position and momentum variables x and p (units-free with K = g, 
[£z,Pfc] = may correspond to the quadrature amplitudes of a single 

electromagnetic mode, i.e., the real and imaginary part of the single mode's 
annihilation operator: a = x + ip. At this stage, it is convenient to consider 
position and momentum eigenstates. We may now replace the Hadamard by a 
Fourier transform, 

•^"l^) position — 1 = I dye ^ |2/)position — \P — ^momentum j (1-13) 
V 71 " J -00 

and the C-NOT gates by appropriate beam splitter operations 3 . The input 
states are taken to be zero-position eigenstates \x = 0). The sequence of beam 
splitter operations Bjk(0) is provided by a network of ideal phase-free beam 
splitters (with typically asymmetric transmittance and reflectivity) acting on 
the position eigenstates as 

Bn(0)\xi, X2) = \x\ sin 6 + X2 cos 6, x\ cos 9 — X2 sin 6) = |a^,x 2 ).(1.14) 



3 A possible continuous- variable generalization of the C-NOT gate is |xi, 0:2) — » \x\,xi + X2), where the 
addition modulo two of the qubit C-NOT, |j/i,j/2) —* \yi,Vl ®V2) with 1/1,3/2 = 0, 1, has been replaced 
by the normal addition. However, for the quantum circuit here, a beam splitter operation as described by 
Eq. (1.14) is a suitable substitute for the generalized C-NOT gate. 
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Now we apply this sequence of beam splitters (making an "iV-splitter"), 

Bn-in(it/4:)B N -2n-i (sin" 1 l/Vi*) x ■ • ■ x B 12 (sin" 1 1/Viv) ,(1.15) 

to a zero-momentum eigenstate \p = 0) oc f dx \x) of mode 1 (the Fourier 
transformed zero-position eigenstate) and N — 1 zero-position eigenstates 
\x = 0) in modes 2 through N. We obtain the entangled iV-mode state 
f dx\x,x, . . . ,x). This state is an eigenstate with total momentum zero and 
all relative positions x\ — Xj = (i, j = 1, 2, . . . , TV). It is clearly an analogue 
to the qubit GHZ state with perfect correlations among the quadratures. How- 
ever, it is an unphysical and unnormalizable state (e.g., for two modes, it cor- 
responds to the maximally entangled, infinitely squeezed two-mode squeezed 
vacuum state with infinite energy). Rather than sending infinitely squeezed 
position eigenstates through the entanglement-generating circuit, we will now 
use finitely squeezed states. 

In the Heisenberg representation, an ideal phase-free beam splitter operation 
acting on two modes with annihilation operators and q is described by 



sin 9 cos 6 \ f Ck 
cos 6 — sin 6 ) I q 



(1.16) 



Let us now define a matrix Bm(0) which is an iV-dimensional identity matrix 
with the entries I^k, Iki, hk, and In replaced by the corresponding entries of 
the above beam splitter matrix. Thus, the matrix for the iV-splitter becomes 

U{N) = B N _ 1N [sm- 1 B N _2N-i (sin- 1 -i= 



x---x J B 12 ^sin- i -^=J . (1.17) 

The entanglement-generating circuit is now applied to N position-squeezed 
vacuum modes. In other words, one momentum-squeezed and N — 1 position- 
squeezed vacuum modes are coupled by an iV-splitter, 

(a; a! 2 ■■■ a' N f=U{N){ 01 a 2 ■■■ a N ) T , (1.18) 
where the input modes are squeezed [13] according to 

ai = coshrid^ + sinhrid^ ^, 

hi = coshr2d l -°' ) — sinh^df"^ , (1.19) 

with z = 2,3, N and vacuum modes labeled by the superscript '(0)'. In 
terms of the input quadratures, we have 

X\ = e L x\ , pi = e L p\ , 

Xi = e-^xf\ p i = e +ra pi 0) , (1.20) 
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for cbj = Xj + ipj (j = 1,2,...,N). The squeezing parameters n and r 2 
determine the degree of squeezing of the momentum-squeezed and the N — 1 
position-squeezed modes, respectively. The correlations between the output 
quadratures are revealed by the arbitrarily small noise in the relative positions 
and the total momentum for sufficiently large squeezing n and ri, 

((4 - x\f) = , {{p\ +p' 2 + .. .+p' N ) 2 ) = iVe-^/4 , (1.21) 

for k ^ I (k, I = 1, 2, N) and a' k = x' k + ip' k . Note that all modes involved 
have zero mean values, thus the variances and the second moments are identical. 

2.3 MEASURING MULTIPARTITE ENTANGLEMENT 

Rather than constructing a circuit for generating entangled states, now our 
task shall be the measurement of multi -party entanglement, i.e., the projection 
onto the basis of maximally entangled multi-party states. For qubits, it is well- 
known that this can be achieved simply by inverting the above entanglement- 
generating circuit (a similar strategy also works for d-level systems [20]). The 
GHZ basis states for N qubits read 

1 

+ (-l) n |l) ® |1 © mi) (g> |1 © m 2 ) <g> • • • (8) |1 mjv-i)), (1.22) 



|*n,mi,m 2 ,...,mjv-i) = — (1°) ® \ m l) ® \ m 2) © • • • © |mjv-l) 



where n, mi, m2, uin-i = 0,1. The projection onto the basis states 
{\^n,m 1 ,m2,...,m N _ 1 }} is accomplished when the output states of the inverted 
circuit (see Fig. 1.1), 

(CNOT JV _ij V CNOT Ar -2JV-i • • • CNOT^^)- 1 

= ^CNOTiaCNOTss ■ ■ ■ CNOTiv-iiv, (1-23) 

are measured in the computational basis. Eventually, {| v I / n,mi,m2,...,mjv-i)} 
are distinguished via the measured output states 

\n) © |mi) © |mi © 171,2) © |m2 © m.3) © • • • © |mAr_2 © mN-i) ■ (1-24) 

Reentering the domain of continuous variables, let us now introduce the 
maximally entangled states 

1 f°° 

\ty(v,ui,U2, ...,un-i)) = —^= / dx e 2lvx \x) © \x — u\) 

J-00 

© \x — u\ — U2) © ••• © \x — ui — U2 — ■ ■ ■ — ujv-i) ■ (1-25) 
Since \x)(x\ = 11 and (x\x') = 5(x — x'), they form a complete, 

/oo 
dvdUl du 2 • • ■dUN-l\ i ^(v,Ui,U2, ...,UN-l)){^(v,Ui,U2, ...,un-i)\ 
-00 

= t® N , (1.26) 
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and orthogonal, 

(*0,Ul,U 2 , ...,ltJV-l)|*(f',«'l,M2, - ViV-l)) 

= <5(v - i/)<5(ui - n / 1 )5(n 2 - u' 2 ) • • • S(un-i - u'n-i)i (I- 27 ) 

set of basis states for N modes. For creating continuous-variable entan- 
glement, we simply replaced the C-NOT gates by appropriate beam split- 
ter operations. Let us employ the same strategy here in order to measure 
continuous-variable entanglement. In other words, a projection onto the 
continuous- variable GHZ basis , iti, «2, un-i))} shall be performed 
by applying an inverse iV-splitter followed by a Fourier transform of mode 
1 and by subsequently measuring the positions of all modes. For an N- 
mode state with modes b\, 62, 6jv, this means that we effectively measure 
Im^ = p' 1: Reb 2 = x 2 ,Re6 3 = x' 3 , ...,Reb N = x' N , with 

( b'i K ■■■ y N ) T =uHN)( h b 2 ••• b N ) T . (1.28) 
For instance, in the three-mode case, the measured observables are 

Pi = -^{Pl + P2 + P3) , 
/2\ 1 

x 2 = y 3^1 - ^/|( x 2 + x 3 ) , 

4 = ^=(z 2 - s 3 ) , (1-29) 

where here bj = xj + ipj. In fact, in a single shot, the quantities vj\f?>, 
y / 2~73(ui + U2/2), and U2/V2 are determined through these measurements, 
and so are all the parameters v = pi +P2+P3, ui = x\ — x 2 , and -u 2 = x 2 — x% 
required to detect the basis state \^>(v, u\, u 2 )) from Eq. (1.25) with N = 3. 
In general, for arbitrary N, the measurements yield p[ = v/y/N and 



/JV-1 / N-2 f N -3 , 



X 



N-3 



x N-2 



x N-l 



4 / 3 / 2 / 1 

- I "Uiv_4 + - I + - I UjV_ 2 + -UN-1 



, , ^AT-3 + 3 I u N-2 + ^ ! 



3 fuiv_2 + 2«JV-] I 



V2 



UN-l, 



(1.30) 
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where v = pi + p 2 H h pn, u± = xi — x 2 , u 2 = x 2 - x 3 ,..., and un~i = 

xn-i—xn. This confirms that the inverse iV-splitter combined with the appro- 
priate homodyne detections (that is tools solely from linear optics) enables in 
principle a complete distinction of the basis states {\^(v,ui,u 2 , ...,un-i)}} 
in Eq. (1 .25). More precisely, the fidelity of the state discrimination can be arbi- 
trarily high for sufficiently good accuracy of the homodyne detectors. We may 
conclude that the requirements of such a "GHZ state analyzer" for continuous 
variables are easily met by current experimental capabilities. This is in contrast 
to the GHZ state analyzer for photonic qubits [capable of discriminating or 
measuring states like those in Eq. (1.22)]. Although arbitrarily high fidelity 
can be approached in principle using linear optics and photon number detec- 
tors, one would need sufficiently many, highly entangled auxiliary photons and 
detectors resolving correspondingly large photon numbers [21, 20]. Neither of 
these requirements is met by current technology. Of course, the C-NOT gates 
of a qubit GHZ state measurement device can in principle be implemented 
via the so-called cross Kerr effect using nonlinear optics. However, on the 
single-photon level, this would require optical nonlinearities of exotic strength. 

In this section, we have shown how measurements onto the maximally 
entangled continuous-variable GHZ basis can be realized using linear optics and 
quadrature detections. These schemes are an extension of the well-known two- 
party case, where the continuous-variable Bell basis [Eq. (1.25) with N = 2] is 
the analogue to the qubit Bell states [Eq. (1.22) with N = 2]. The continuous- 
variable and the qubit Bell states form those measurement bases that were used 
in the quantum teleportation experiments [22] and [23, 24], respectively. The 
extension of measurements onto the maximally entangled basis to more than 
two parties and their potential optical realization in the continuous-variable 
realm might be relevant to multi-party quantum communication protocols such 
as the multi-party generalization of entanglement swapping [25]. However, the 
entanglement resources in a continuous-variable protocol, namely the entangled 
continuous-variable states that are producible with squeezed light and beam 
splitters, exhibit only imperfect entanglement due to the finite degree of the 
squeezing. When can we actually be sure that they are multi-party entangled at 
all? In the next section, we will address this question and discuss criteria for the 
theoretical and the experimental verification of multipartite continuous-variable 
entanglement. 

2.4 SUFFICIENT INSEPARABILITY CRITERIA 

For continuous-variable two-party states, an inseparability criterion can be 
derived that does not rely on the partial transpose. It is based on the variances 
of quadrature combinations such as x\ — x 2 and p\ + p 2 , motivated by the fact 
that the maximally entangled bipartite state J dx \x,x) is a (zero-)eigenstate 
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of these two combinations [26]. Similarly, in a continuous- variable Bell mea- 
surement, the quadrature combinations x\ — x 2 and p\ + p 2 are the relevant 
observables to be detected. Hence, for two modes, applying the variance- 
based inseparability criterion and measuring in the maximally entangled basis 
can both be accomplished by equal means, namely a single beam splitter and 
two homodyne detectors. In other words, the effectively inverse circuit for 
the generation of bipartite entanglement provides the recipe for both measur- 
ing maximum entanglement and verifying nonmaximum entanglement. When 
looking for multi-party inseparability criteria for arbitrarily many modes, it 
seems to be natural to pursue a similar strategy. We are therefore aiming at a 
criterion which is based on the variances of those quadrature combinations that 
are the measured observables in a continuous-variable GHZ measurement. 
Let us consider three modes. According to Eq. (1.29), we define the operators 

u = -j=(x 2 -x 3 ), 

d = vf £l ~7i ( * 2 + * 3) ' 

w = —7={pi +P2 +P3) x \/2 , (1.31) 
v3 

where we added a factor of y/2 in w compared to the first line of Eq. (1.29). 
Let us further assume that the three-party state p is fully separable and can be 
written as a mixture of tripartite product states, 



Pil ® Pi2 ® Pi3 ■ (1-32) 



Using this state, we can calculate the total variance of the operators in Eq. (1.31), 



((Auf) p + ((Av) 2 ) p + ((Aw) 2 ) p 



= £ P t ((u% + (v% + (w%) - (u) 2 p - {v) 2 p - (w)l 

i 

= E \ «*i>< + < £ 2>» + (4h + (Pi>< + (pl)i + (Pih) 

i 

-^Pi o (ixi)i{x 2 )i + (Xl)i(x 3 )i + (X 2 )i{x3)i 

i 

- 2(pi)i(p2)i - 2(pi) i (p 3 ) i - 2(p 2 ) 4 (p 3 ) l ) - (u) 2 p - (v) 2 p - (w) 



r.\2 
P 
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Y,Pi\ (<(Axi) 2 >, + ((Ax 2 ) 2 ) t + ((Ax 3 )% 



i 



+ ((Ap 1 ) 2 ) t + ((Ap 2 ) 2 ) t + ((Ap 3 ) 2 ) l ) 

+E P * < fi >? - (e p * <*>«) + E p * - fe p * 




2 



i \ i / i \ i 



+E p *<< B >. ? - E^w* 



(1.33) 



i \ i 



where (•••)« means the average in the product state pn <g) fa <8> fa- Similar 
to the derivation in Ref. [26], we can apply the Cauchy-Schwarz inequality 
Y^i Pi{u)i > (Yli Pi\(u)i\) 2 , and see that the last two lines in Eq. (1.33) are 
bounded below by zero. Also taking into account the sum uncertainty relation 
{{Axj) 2 )i + {(&Pj) 2 )i > \[xj,Pj]\ = 1/2 (j = 1, 2, 3), we find that the total 
variance itself is bounded below by 1 (using J2i Pi = !)• ^ n y tota l variance 
smaller than this boundary of 1 would imply that the quantum state concerned is 
not fully separable as in Eq. (1.32). But would this also imply that the quantum 
state is genuinely tripartite entangled in the sense that none of the parties can 
be separated from the others (as, for example, in the pure qubit states |GHZ) 
and | W))? This is obviously not the case and a total variance below 1 does not 
rule out the possibility of partial separability. The quantum state might still 
not be a genuine tripartite entangled state, since it might be written in one or 
more of the following forms 4 [18]: 



P = E Pi P il2 &3, P = E P 'i P il3 P = E P i P i23 ® Pa- (1-34) 



Thus, in general, a violation of ((An) 2 ) p + ((Av) 2 ) p + ((Aw) 2 } p > 1 does not 
necessarily witness genuine tripartite entanglement (a counterexample will be 
given below). However, it does witness genuine tripartite entanglement when 
the quantum state in question is pure and totally symmetric with respect to all 
three subsystems [28]. In that case, a possible separation of any individual 
subsystem, 



implies full separability, p = p\ (g> p2 pz- Hence a total variance below 1 
negates the possibility of any form of separability in this case. 



P = P12 ® P3, P = P13 ® P2, P = P23 <8> Pi 



(1.35) 



4 A full classification of tripartite Gaussian states is given in Ref. [27] in analogy to that for qubits from 
Ref. [18]. In addition, necessary and sufficient three-mode inseparability criteria for Gaussian states are 
proposed in Ref. [27]. 
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By extending the quadrature combinations in Eq. (1.31) from 3 to N parties 
(corresponding to the output modes of an inverse iV-splitter) and performing 
a similar calculation as for N = 3 with an additional factor of y/N — 1 in 
the total momentum operator w, we find that any iV-mode state with modes 
b\, b 2 , •■•) t>N which is fully separable, p = ^ Pj pn ® p i2 (8> • • • <8> Pin, obeys 
the inequality 

((Apl) 2 ), + Ef =^f 2) ^ > ± • (1-36) 

Here, p[ = Imb^x^ = Reb 2 , ...,x f N = Reb' N are the corresponding out- 
put quadratures of the inverse iV-splitter applied to the modes b\, b 2 , 1>n 
[Eq. (1.28)]. Alternatively, one can also derive the following necessary condi- 
tion for full separability [28], 



2(N-1) >^ 2 



+ ((AP)^ > - . (1.37) 



In this inequality, Xij = Xi — Xj and P = J2iLi Pi we tne relative positions 
and the total momentum of the relevant state with modes bj = Xj + ipj. 

The choice of the operators in the inequalities Eq. (1.36) and Eq. (1.37) relies 
upon the fact that the quantum fluctuations of these observables simultaneously 
vanish for maximum GHZ entanglement. This must be in agreement with their 
commutation relations, and indeed, we have 

[Xtj , P] = [xi - Xj , pi + pj] = , (1.38) 

for any N and i, j. Similarly, the output quadratures of the inverse iV-splitter 
yield, for instance, for N = 3, 

[Pi +P2+ fa, x 2 - x 3 ] = 0, [pi +P2+ P3, 2*i - (x 2 + x 3 )\ = . (1.39) 

Both criteria in Eq. (1.36) and Eq. (1.37) represent necessary conditions for 
full separability, though they are not entirely equivalent [i.e., there are partially 
inseparable states that violate Eq. (1.37), but satisfy Eq. (1.36), see below]. 
Moreover, the criterion in Eq. (1.37) contains in some sense redundant observ- 
ables. As we know from the previous section, N observables suffice to measure 
an iV-party GHZ entangled state. These N observables are suitably chosen 
quadratures of the N output modes of an inverse iV-splitter. Their detection 
simultaneously determines the total momentum and the N — 1 relative positions 
£1 — x 2 , x 2 — £3,..., and £jv-i — £jv [Eq. (1.30)] 5 . However, in Eq. (1.37), 



5 The variances of the N — 1 relative positions x\ — £2, X2— £3 andxjv-i — xn are also available via the 

variances of the output quadratures of the inverse TV-splitter. First, the variance of£^ = (£ jv - 1 — i jv ) 
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there are 1 + [N(N — l)]/2 different operators. Nevertheless, for two parties 
and modes, the conditions in Eq. (1.36) and Eq. (1.37) coincide and correspond 
to the necessary separability condition for arbitrary bipartite states given in 
Ref. [26]. 

In summary, we have shown in this section that the circuit for measur- 
ing GHZ entanglement also provides a sufficient inseparability criterion for 
arbitrary multi-party continuous-variable states (pure or mixed, Gaussian or 
non-Gaussian) of arbitrarily many parties. This criterion is experimentally ac- 
cessible via linear optics and homodyne detections. The disadvantage of not 
being a necessary inseparability condition (not even for Gaussian states, see 
below) might be unsatisfactory from a theoretical point of view, but would not 
be an obstacle to experimental inseparability proofs. A more serious drawback, 
in particular when considering an experimental verification of multi-party en- 
tanglement, is the fact that without additional assumptions, for arbitrary states, 
the criteria presented in this section are in general not sufficient for genuine 
multi-party inseparability. They verify only partial inseparability. In the next 
section, we will give a simple example for this. 



Partial multipartite entanglement. Let us investigate how the following 
pure three-mode state described by the Heisenberg operators 



behaves with respect to the multi-party inseparability criteria. Modes 1 and 2 
are in a two-mode squeezed vacuum state [Eq. (1.18) and Eq. (1.20) for N = 2 
and r = r± = r^\, and mode 3 is in the vacuum state. This three-party state 
is obviously only partially 6 entangled (it is an example for the second class 
of the five classes of three-mode Gaussian states in Ref. [27]). Applying an 
inverse "tritter" (three-splitter) to these modes, calculating the relevant output 



corresponding to the last line in Eq. (1.30) is directly measurable. In addition, by converting the measured 
photocurrent into a light amplitude and "displacing" (feed-forward) &' N _ 1 according to x' N _ 1 —* x'^_ l = 

x' N _ 1 — -^ X ' N = y §(£jv-2 — Xjv_i), one can directly measure the variance of x N _ 2 — x N-i etc - 
Similarly, one would also employ this feed-forward technique in a multi-party quantum communication 
protocol that relies on the N classical results of an Af-mode GHZ state measurement. 
6 note that we use the term "partially entangled" here for states which are not genuinely multi-party insepa- 
rable. In the literature, sometimes "partial entanglement" is also referred to as nonmaximum entanglement 
of two or more parties (in the sense that for two parties the Schmidt coefficients are not all equal). As 
discussed later, also the genuinely multi-party entangled continuous-variable states are only nonmaximally 
entangled due to the finite degree of the squeezing. 



2.5 



MULTI-PARTY ENTANGLED STATES 




(1.40) 
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Figure 1.2 Application of the necessary conditions for full three-party separability. On the left: 
the inequality Eq. (1.42) with the boundary 3/2 as a function of the squeezing r; for < r < 1, 
the inequality is violated nearly everywhere. On the right: the inequality Eq. (1.41) with the 
boundary 1/2 as a function of the squeezing r; for < r < 1, the inequality is statisfied nearly 
as much as it is violated. 



variances, and inserting them into Eq. (1.36) yields 

l(l e +^ + e -A + l>i (i.4i) 

4V3 / 6 - 2 ' 

as a necessary condition for full separability. We find that equality holds for 
r = which doesn't tell us anything, though we know, of course, that the state is 
fully separable in this case and must obey Eq. (1.36). For some finite squeezing 
r, the inequality Eq. (1.41) is illustrated in Fig. 1.2. Similarly, application of 
the criterion in Eq. (1.37) to the above state leads to 

1 3 

- (3e~ 2r + cosh 2r + 2) > - . (1.42) 

Again, equality holds for r = 0. In Fig. 1.2, also this condition is depicted for 
some finite squeezing r. 

The comparison between the two conditions in this example demonstrates 
that they are not equivalent. For some squeezing, the partially entangled three- 
mode state violates Eq. (1.37) while satisfying Eq. (1.36). Moreover, both 
conditions can apparently be violated by an only partially entangled state. 
Hence, both their violation, though ruling out full separability, does not imply 
the presence of genuine multi-party entanglement. Another observation is that 
both conditions are satisfied for sufficiently large squeezing when the partial 
entanglement is sufficiently good. This confirms that the two conditions are 
necessary for full separability, but not sufficient, not even for a Gaussian state 
like that in our example. In fact, also the bipartite separability condition 
of Ref. [26] is both necessary and sufficient only for Gaussian states in a 
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certain standard form (where any Gaussian state can be transformed into this 
standard form via local operations). The partially entangled three-mode state 
here lacks the symmetry that is required for a state to violate the separability 
conditions whenever it contains some entanglement. We will now turn to a 
family of multipartite entangled states which are totally symmetric with respect 
to all their parties, which do always violate both conditions for full multi-party 
separability, and which are indeed genuinely multi-party entangled. 

Genuine multipartite entanglement. Let us consider the output states of the 
entanglement-generating circuit in section 2.2. There, we applied an iV-splitter 
to one momentum-squeezed (r{) and N — 1 position-squeezed (r 2 ) vacuum 
modes to obtain the modes a^a^, ■ ■■,a' N [Eq. (1.18)]. Now one can easily 
see that the first multi-party separability condition is violated for any nonzero 
squeezing r\ > or r 2 > 0, because application of an inverse iV-splitter means 



(a? a 2 ' ••• a" N f = U\N)(a[ a' 2 ••• a' N f (1.43) 

= U\N)U(N) ( ai a 2 ••• a N ) T , 



with ai, a,2, ajv fr° m Eq. (1.19). Since W(N)U(N) = I (identity matrix), 
the squeezed quadratures of Eq. (1.20) can be directly inserted into Eq. (1.36) 
yielding a violation of (e~ 2ri + e _2r2 )/4 > 1/2 for any n > or r 2 > 0. 
Alternatively, using Eq. (1.21), the criterion in Eq. (1.37) becomes 



This condition is also violated for any r\ > or r 2 > 0. Due to their purity 
and total symmetry we conclude that the members of the family of states which 
emerge from the iV-splitter circuit are genuinely multi-party entangled for any 
r\ > or r 2 > 0. This applies in particular to the case where r\ > and 
T2 = 0, i.e., when only one squeezed light mode is required for the creation of 
genuine multipartite entanglement. 

Independent of the inequalities Eq. (1.36) and Eq. (1.37), there are also other 
ways to see that these particular states are genuinely multi-party entangled. One 
simply has to find some form of entanglement in these states. For example, by 
tracing out modes 2 through N of the pure iV-mode state given in Eq. (1.18), 
one finds that the remaining one-mode state is mixed, provided r\ > or r 2 > 
[28]. Thus, the pure iV-mode state is somehow entangled and hence genuinely 
multi-party entangled due to its complete symmetry. Note that in order to infer 
even only partial entanglement via tracing out parties, the iV-mode state here 
has to be pure. In contrast, the multi-party inseparability criteria of section 2.4 
may verify partial inseparability for any iV-mode state. 

From a conceptual point of view, it is very illuminating to analyze which 
states of the above family of iV-mode states can be transformed into each other 




(1.44) 
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via local squeezing operations [29]. For example, by applying local squeezers 
with squeezing si and s 2 to the two modes of the bipartite state generated with 
only one squeezer [Eq. (1.20) for N = 2 with r 2 = 0], we obtain 



x'{ 


= e- Sl x[ 


= (e 


p'l 


= e+Vi 


= (e 


x 2 


= e~ S2 x' 2 


= (e 


p'i 


= e +S2 p' 2 


= (e 



+-i-^ )+ e +«i4°))/v^, 

- e -«4 0) )/V^, 
+ S2 -r lj5 (o) _ e + S 2^(0))/^2 . (1.45) 

With the choice of si = s 2 = r\/2 = r, the state in Eq. (1.45) is identical to 
a two-mode squeezed state built from two equally squeezed states [Eq. (1.20) 
for N = 2 with n = r, r 2 = r]. The latter and the state produced with only 
one squeezer [Eq. (1.20) with r\ = 2r and r 2 = 0] are equivalent under local 
squeezing operations. This means that Alice and Bob sharing the state produced 
with one squeezer r\ = 2r have access to the same amount of entanglement as 
in the "canonical" two-mode squeezed state with squeezing r = r\j% -E V .N. = 
[cosh(r 1 /2)] 2 log[cosh(ri/2)] 2 - [sinh(ri/2)] 2 log[sinh(n/2)] 2 [30]. For a 
given amount of entanglement, however, the canonical two-mode squeezed 
vacuum state has the least mean photon number. Conversely, for a given mean 
energy, the canonical two-mode squeezed vacuum state contains the maximum 
amount of entanglement possible. 

Similar arguments apply to the states of more than two modes. From the 
family of iV-mode states, the state with the least mean photon number is 
determined by the relation 



e ±2ri = (iv- l)sinh2r 2 



1 + (N - l) 2 sinh 2 2r 2 ± 1 



(1.46) 



This relation is obtained by requiring each mode of the iV-mode states to be 
symmetric or "unbiased" in the x and p variances [29]. Only for iV = 2, we 
obtain n = r 2 . Otherwise, the first squeezer with n and the N — 1 remaining 
squeezers with r 2 have different squeezing. In the limit of large squeezing, we 
may use sinh2r 2 rs e +2r2 /2 and approximate e +2ri of Eq. (1.46) by 

e +2ri ~ (N — l)e +2r2 . (1.47) 

We see that in order to produce the minimum-energy iV-mode state, the single 
ri-squeezer is, in terms of the squeezing factor, N — 1 times as much squeezed 
as each r 2 -squeezer. However, also in this general iV-mode case, the other 
iV-mode states of the family can be converted into the minimum-energy states 
via local squeezing operations. This applies in particular to the A^-mode states 
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produced with just a single squeezer and to those built from N equally squeezed 
states. As a result, due to the equivalence under local entanglement-preserving 
operations, with a single sufficiently squeezed state and beam splitters, arbi- 
trarily many genuinely multi -party entangled modes can be created just as well 
as with N squeezers and beam splitters. 

In contrast to the three-mode state given by Eq. (1.40), the output states of 
the iV-splitter are totally symmetric under interchange of parties. This becomes 
more transparent when we look at the states in the Wigner representation. For 
simplicity, let us assume r = n = ri- The position-squeezed input states of 
the iV-splitter circuit, for instance, have the Wigner function 

W(x,p) = - exp(-2e +2r x 2 - 2e~ 2r p 2 ). (1.48) 

7T 

Through the linear iV-splitter operation, the total input Wigner function to the 
iV-splitter (one momentum-squeezed and N — 1 position-squeezed vacuum 
modes), 

f2\ N 

Win(x,p) = I -J exp(-2e- 2r x 1 2 - 2e +2r Pl 2 ) (1.49) 

exp(-2e +2r x 2 2 - 2 e - 2r p 2 2 ) exp(-2e+ 2r x 3 2 - 2e~ 2r p 3 2 ) 
x • • • x exp(-2e +2r X7v 2 - 2eT 2T p N 2 ), 

Pi? 



x 3 ) 2 



(1.50) 

Here we have used x = (xi,X2, —,xn) and p = (jpi,P2, ■ ■■,Pn)- The pure- 
state Wigner function W ou t(x, p) is always positive, symmetric among the N 
modes, and becomes peaked at Xj — x 3 ■ = = 1,2, ...,N) and p\ + 
P2 + • • • + Pn = for large squeezing r. For N = 2, it exactly equals the 
well-known two-mode squeezed vacuum state Wigner function [13], which is 
proportional to S(xi — X2)S(pi + P2) in the limit of infinite squeezing. As 
discussed previously, the state W / out (x, p) is genuinely iV-partite entangled for 
any squeezing r > 0. The quantum nature of the cross correlations x; L x 3 - and 
PiPj appearing in W ut(x, p) for any r > is also confirmed by the purity of 



is transformed into the output Wigner function 
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this state. This purity is guaranteed, since beam splitters turn pure states into 
pure states (it can also be checked via the correlation matrix of the Gaussian 
state W ut(x,p) [28]). 

A nice example for a multi-party entangled state which is not a member of 
the above family of states and hence not totally symmetric with respect to all 
its modes is the (M + l)-mode state described by the Wigner function 



/ 2 \M+i r / cog M+i 

Wmqc(x,p) = i^-j exp | - 2e~ 2ri ( sin0 o xi + —j=- ^ 



2 



Xi 



/ . a M+1 \ 2 

-2e+ 2 ^ cos^ox!-^ Y xA 
V ^ U J 

i M+1 1 

(1.51) 



where x = (x 1 ,x 2 , ...,x M +i), P = (pi,P2, -,Pm+i), and 

- 1 < singp < J ^ - , (1.52) 

-2n _ sm #o - cos 6*o _2r 2 _ V^Mcos^o — sin 6> 
\/M sin 6<o + cos 6q ' cos 6q + sin 6q 

The significance of this (M + l)-mode state is that it represents a kind of 
multiuser quantum channel ("MQC") enabling optimal 1 — > M "telecloning" 
of arbitrary coherent states from one sender to M receivers [31]. Though 
not completely symmetric with respect to all M + 1 modes (but to modes 2 
through M + 1), it is a pure Gaussian state which is indeed genuinely multi- 
party entangled. This can be seen, because none of the modes can be factored 
out of the total Wigner function. Despite its "asymmetry", this state is not 
only partially multi-party entangled as is the asymmetric pure three-mode state 
given by Eq. (1.40). Of course, the bipartite entanglement between mode 1 on 
one side and modes 2 through M + 1 on the other side is the most important 
property of VFmqc( x ) p) in order to be useful for 1 — > M telecloning [31]. 
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The generation of the state Wmqc( x j p) is ver y similar to that of the above 
family of multi-party entangled states produced with an iV-splitter: first make 
a bipartite entangled state by combining two squeezed vacua, one squeezed in 
p with r\ and the other one squeezed in x with r 2 , at a phase-free beam splitter 
with reflectivity/transmittance parameter 9 = 9q. Then keep one half (the 
mode 1) and send the other half together with M — 1 vacuum modes through an 
M-splitter. The annihilation operators of the initial modes aj before the beam 
splitters, j = 1, 2, M + 1, are then given by 

a\ = coshria^ + sinhria^ ^, 
d 2 = coshr^dj, '' — sinhr^dj, ^, 
a, = af ] , (1.53) 

where i = 3,4, ...,M + 1. 

By using the ideal phase-free beam splitter operation from Eq. (1.16), with 
Bki(6) this time representing an (M + 1) -dimensional identity matrix with the 
entries I^k, IkU hk, and In replaced by the corresponding entries of the beam 
splitter matrix in Eq. (1.16), the MQC -generating circuit can be written as 



( h b 2 ■■■ b M +i ) = Umqc(M + 1) ( ai a 2 ■■■ a M +i ) • 



(1.54) 



U M qc(M + 1) = Bmm+i ^sin 1 B M -im ^sin 1 -^=\ 

x-x^3 4 (sin-^= T ) 523 (sin--^) 



x 



Bi 2 (0 o ) • (1-55) 



The first beam splitter, acting on modes di and d 2 , has reflectivity/transmittance 
parameter 9 = Oq. The remaining beam splitters represent an M-splitter. In 
Eq. (1 .54), the output modes bj correspond to the M+ 1 modes of the MQC state 
described by Wmqc in Eq. (1.51). Let us now return to the totally symmetric 
multipartite entangled states given by Eq. (1.18) and explore some of their 
properties. For simplicity, we will thereby focus on those states emerging from 
the iV-splitter circuit which are created with input states equally squeezed in 
momentum and position, r = r\ = r 2 . 

Nonlocality and other properties. In this paragraph, we will discuss some of 
the properties of the state W ou t( x ) p) m Eq. (1.50). This will further illustrate 
the character of W ou t ( x > p) as a nonmaximally entangled multi-party state. 
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One such property is that this state, despite having an always positive Wigner 
function, violates iV-party Bell-type [2] (or Mermin-type [32]) inequalities 
imposed by local realism for any squeezing r > [33]. The observables 
producing these violations are displaced photon-number parities rather than 
continuous variables such as x and p [34]. Like for the qubit states [32], 
the violations increase as the number of parties N grows. However, this 
increase becomes steadily smaller for larger N, as opposed to the exponential 
increase for the maximally entangled qubit GHZ states [32]. This discrepancy 
may be explained by the fact that the violations are exposed only for finite 
squeezing where the state W out (x, p) is a nonmaximally entangled multi-party 
state [33]. Note that, in general, the violations of iV-party inequalities imposed 
by local realism do not necessarily imply the presence of genuine multipartite 
entanglement. However, for the pure and symmetric states H / out (x, p), once 
again, proving some kind of entanglement means proving genuine multipartite 
entanglement. 

In order to prove the nonlocality exhibited by the state W (x, p) = W ou t ( x > p) , 
let us now use the fact that the Wigner function is proportional to the quantum 
expectation value of a displaced parity operator [35, 34]: 

W(a) = 0) (n(a)) = 0) n(a) , (1.56) 

where a = x + ip = (a>\, c*2, aw) and U(ct) is the quantum expectation 
value of the operator 

N N 

fi(a) = (g) Ui(ai) = (g) A (<*i) • (1-57) 
i=\ i=i 

The operator Di(cti) is the displacement operator, 

D(a) = exp(cra^ — a* a) , (1.58) 
acting on mode i. Thus, 11(a) is a product of displaced parity operators given 

by 

ftite) = nS +) (oi) - nS _) (ai) , (1.59) 
with the projection operators 

oo 

D i (a i )J2\^)(2k\D\(a i ), (1.60) 

k=0 

oo 

k=0 



ni _) («i) = 
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corresponding to the measurement of an even (parity +1) or an odd (parity —1) 
number of photons in mode i. This means that each mode is now characterized 
by a dichotomic variable similar to the spin of a spin- 1/2 particle or the single- 
photon polarization. Different spin or polarizer orientations from the original 
qubit based Bell inequality are replaced by different displacements in phase 
space. This set of two-valued measurements for each setting is just what we 
need for the nonlocality test. 

In the case of iV-particle systems, such a nonlocality test is possible using the 
iV-particle generalization of the two-particle Bell-CHSH inequality [16]. This 
inequality is based on the following recursively defined linear combination of 
joint measurement results (in this paragraph, the symbol B does not refer to a 
beam splitter operation), 

B N = -[a(a N ) +a(a' N )]B N -i 

+\[a(a N ) - a(a' N )]B' N ^ = ±2 , (1.62) 

where <r(ajv) = ±1 and <r(a' N ) = ±1 describe two possible outcomes for two 
possible measurement settings (denoted by a at and a' N ) of measurements on 
the iVth particle. Note, the expressions B' N wee equivalent to B^ but with all 
the cii and swapped. Provided that Bn-i = ±2 and B' N l = ±2, Equation 
(1.62) is trivially true for a single run of measurements where <j(o,n) is either 
+ 1 or —1 and similarly for a(a' N ). Induction proves Eq. (1.62) for any N 
when we take 

B 2 = [o-(ai) + a(a' l )]a{a 2 ) 

+[a{a 1 )-a{a' 1 )}a{a! 2 ) = ±2. (1.63) 

Within the framework of local realistic theories with hidden variables A = 
(Ai, A2, Ajv) and the normalized probability distribution P(X), we obtain 
an inequality for the average value of Bn = Bn(X), 



/ 



d\id\ 2 ...d\ N P(X)B N (X) 



< 2 . (1.64) 



By the linearity of averaging, this is a sum of means of products of the cr(aj) 
and cr(a'j). For example, if N = 2, we obtain the CHSH inequality 

\C{ ai ,a 2 ) + C(oi, a! 2 ) + C(a[,a 2 ) - C(a\, a' 2 )\ < 2, (1.65) 

with the correlation functions 

C(ai,a 2 ) = J d\ 1 dX 2 P(X 1 ,X 2 )a(a 1 ,Xi)a(a 2 ,X 2 ) . (1.66) 
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Following Bell [2], an always positive Wigner function can serve as the hidden- 
variable probability distribution with respect to measurements corresponding 
to any linear combination of x and p. In this sense, the finitely squeezed two- 
mode squeezed state Wigner function could prevent the CHSH inequality from 
being violated when restricted to such measurements: W(xi,pi,X2,P2) = 
P(X\, ^2)- The same applies to the Wigner function in Eq. (1.50): W(x, p) = 
P(X) could be used to construct correlation functions 

C(a) = J dX 1 dX 2 ...dX N P(\) 

x <r(ai, Xi)a(a 2 , A 2 ) ■ ■ ■ a(a N , X N ), (1.67) 

where a = (a\, a 2 , a/v). However, for parity measurements on each mode 
with possible results ±1 for each differing displacement, this would require 
unbounded <5-functions for the local objective quantities a (en, Aj) [34], as in 
this case we have 

C(a) = n(a) = (71-/2)^(0) . (1.68) 

This relation directly relates the correlation function to the Wigner function 
and is indeed crucial for the nonlocality proof of the continuous- variable states 
in Eq. (1.50). 

Let us begin by analyzing the nonlocal correlations exhibited by the entan- 
gled two-party state. For this state, the two-mode squeezed state in Eq. (1.50) 
with N = 2, we may investigate the combination [34] 

B 2 = n(0, 0) + n(0, (3) + U(a, 0) - U(a, f3) , (1.69) 

which according to Eq. (1.65) satisfies \B 2 \ < 2 for local realistic theories. 
Here, we have chosen the displacement settings a\ = a 2 = and a[ = a, 

ol 2 = 13. 

Writing the states in Eq. (1.50) as 
f N 

11(a) = exp< - 2 cosh 2r^ | ai \ 2 (1.70) 



+ sinh 2r 



i=i 

TV N 



=1 



for N = 2 and a = (3 = i\fj with the real displacement parameter J > 1 , 
we obtain 5 2 = 1+2 exp(-2J" cosh 2r)-exp(-4 l 7e +2r ). In the limit of large 



'This choice of two equal settings leads to the same result as that of Banaszek and Wodkiewicz [34] who 
used opposite signs: a = \fj and = —\f~3. 
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r (so cosh 2r m e +2r /2) and small J", £> 2 is maximized for Je +2r = (In 2)/3, 
yielding S™ ax « 2.19 [34], which isaclear violation of the inequality |,B 2 | < 2. 
Smaller violations also occur for smaller squeezing and larger J. Indeed, for 
any nonzero squeezing, some violation takes place [33]. 

We will now consider more than two parties. Let us first examine the three- 
mode state by setting N = 3 in Eq. (1.50). According to the inequality of the 
correlation functions derived from Eq. (1.62)-(1.64), we have 

\C(ai, a 2 ,a 3 ) + C(ai,a' 2 ,a 3 ) + C(a[,a 2 ,a 3 ) - C(ai, a 2 , a' 3 )| < 2. (1.71) 



Thus, for the combination 



B 3 = 11(0, 0, 7) + 11(0, 0, 0) + n(a, 0, 0) - 1%, P, 7), (1.72) 

a contradiction to local realism is demonstrated by |^3| > 2. The corresponding 
settings here are ot\ = a 2 = a 3 = and a[ = a, a' 2 = (3, a f 3 = 7. With the 
choice a = y/Je** 1 , (3 = yfje i<t>2 , and 7 = y/je l<t> ' A , we obtain 

3 2 
B 3 = ^ exp (- 2.7 cosh 2r - -J" sinh 2r cos 2^) (1.73) 



i=l 



— exp < 



1 

- 6i7 cosh 2r — - J sinh 2r [cos 20j — 4 cos (<^ + (f>j ) 



Apparently, because of the symmetry of the entangled three-mode state, 
equal phases should also be chosen in order to maximize B 3 . The best 
choice is 0! = </> 2 = 4> 3 = it/ 2, which ensures that the positive terms in 
Eq. (1.73) become maximal and the contribution of the negative term minimal. 
Therefore, we again use equal settings a = (3 = 7 = i\fj and obtain 

B 3 = 3exp(-2Jcosh2r + 2^sinh2r/3) - exp(-6^e +2r ) . (1.74) 

The violations of | B 3 \ < 2 that occur with this result are similar to the violations 
of \B 2 \ < 2 obtained for the two-mode state, but the N = 3 violations are even 
more significant than the N = 2 violations [33]. In the limit of large r (and 
small J), we may use cosh2r w sinh2r ps e +2r /2 in Eq. (1.74). Then B 3 
is maximized for J e +2r = 3(ln3)/16: Bf^ w 2.32. This optimal choice 
requires smaller displacements J than those of the N = 2 case for the same 
squeezing. 

Let us now investigate the cases N = 4 and N = 5. From Eq. (1.62)- 
(1.64) with N = 4, the following inequality for the correlation functions can 
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be derived: 

^|C(ai,o 2 ,o 3 ,a 4 ) + (7(ai, a 2 ,o 3 ,a 4 ) + C(ai, a' 2 , a 3 , a 4 ) 

+C(ai,a 2 ,a 3 ,a4) + C(ai, a 2 , a' 3 , a 4 ) + C(ai, a 2 , a 3 , a 4 ) 
+(7(0^,02,03,04) + C(ai,a 2 ,a 3 ,a 4 ) + C(a' l5 a2, a' 3 , 04) 
+(7(ai, a' 2 , 03, a 4 ) — (7(0^, a' 2 , a' 3 , 04) — (7(6^, a 2 , 03, a' 4 ) 
-(7(^,02, 03,04) - C(ai, a 2 ,a 3 ,a 4 ) - C(oi, o 2 , a 3 , a 4 ) 
-C(ai, 02,03,04)! < 2 . (1.75) 

It is symmetric among all four parties as any inequality derived from Eq. (1 .62)- 
(1.64) is symmetric among all parties. For the settings ot\ = 02 = 03 = 04 = 
and a[ = a, a 2 = /3, a' 3 = 7, a 4 = (5, complying with local realism means 
|B 4 | < 2 where 

£ 4 = ^n(o,o,o,o) + n(o,o, 7 ,o)+n(o,/3,o,o) 

+n(a, 0, 0, 0) + n(o, o, 7 , <5) + n(o, p, o, <5) 
+n(a, 0, 0, 5) + n(o, p, 7, 0) + II(a, o, 7, 0) 

+II(a, (3, 0, 0) - n(a, (3, 7, 0) - n(a, (3, 0, 5) 
-n(a, 0, 7, 5) - n(o, /3, 7, <5) - n(o, 0, 0, 0) 
-n(a,/?, 7 ,«5)] . (1.76) 

Similarly, for iV = 5 one finds 

B 5 = i [n(o, 0, 0, 5, e) + n(o, 0, 7, 0, e) + n(o, /?, 0, 0, e) 

+n(a, 0, 0, 0, e) + n(o, 0, 7, 5, 0) + n(0, f3, 0, (5, 0) 
+n(a, 0, 0, 5, 0) + n(o, p, 7, 0, 0) + II(a, 0, 7, 0, 0) 
+n(a, /3, 0, 0, 0) - II(a, /3, 7, 5, 0) - n(a, 7, 0, e) 
-n(a, /3, 0, 5, e) - U{a, 0, 7, 5, e) - 11(0, P, 7, 5, e) 
-n(0,0,0,0,0)] , (1.77) 

which has to statisfy \B$\ < 2 and contains the same settings as for N = 4, but 
in addition we have chosen 05 = and a' 5 = e. 

We can now use the entangled states of Eq. (1.70) with N = 4 and N = 5 
and apply the inequalities to them. For the same reason as for N = 3 (symmetry 
among all modes in the states and in the inequalities), the choice a = (3 = 7 = 
5 = e = %\fj appears to be optimal (maximizes positive terms and minimizes 
negative contributions). 
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With this choice, we obtain 

B 4 = 2exp(-2 l 7cosh2r + Jsmh2r) 

-2 exp(-6.7 cosh 2r - 3 J sinh 2r) 

+3exp(-47cosh2r) - ^ exp(-8.7e +2r ) - - , 

B 5 = 5exp(-4 l 7cosh2r + 4 l 7sinh2r/5) 

5 1 
-- exp(-8.7 cosh 2r - 247 sinh 2r/5) - - . (1.78) 

Apparently, the maximum violation of \Bn\ < 2 (for our particular choice 
of settings) grows with increasing number of parties N [33]. The asymptotic 
analysis (large r and small J) yields, for instance, for N = 5: B" iax f» 
2.48 with Je +2r = 5(ln2)/24. For a given amount of squeezing, smaller 
displacements J than those for N < 4 (at the same squeezing) are needed to 
approach this maximum violation. Another interesting observation is that in 
all four cases (N = 2, 3, 4, 5), violations occur for any nonzero squeezing [33]. 
This implies the presence of iV-partite entanglement for any nonzero squeezing. 
Moreover, also for modest finite squeezing, the size of the violations (at optimal 
displacement J) grows with increasing N [33]. 

Larger numbers of parties N > 5 were also considered in Ref. [33]. The 
degree of nonlocality of the continuous-variable states, if represented by the 
maximum violation of the corresponding Bell-type inequalities, seems to grow 
with an increasing number of parties. This growth, however, decelerates for 
larger numbers of parties. Thus, the 'evolution' of the continuous-variable 
states' nonlocality with an increasing number of parties and the corresponding 
'evolution' of nonlocality for the qubit GHZ states are qualitatively similar but 
quantitatively different with an exponential increase for the qubits. The reason 
for this may be that the qubit GHZ states are maximally entangled, whereas the 
continuous- variable states are nonmaximally entangled for any finite squeezing. 
Similarly, the N -party version of the nonmaximally entangled qubit state |W) 
yields a non-exponential increase of the maximum violations (by employing, 
for example, an analysis analogous to that here [36]). Note that the observation 
of the nonlocality of the continuous-variable states here requires small but 
nonzero displacements J oc e~ 2r , which is not achievable when the singular 
maximally entangled states for infinite squeezing are considered. 

Finally, the "unbiased" minimum-energy states of the family of entangled 
iV-mode states might yield larger violations. These states are not produced with 
N equal squeezers (as those states whose nonlocality we have analyzed here), 
but with one ri-squeezer and N — 1 ^-squeezers related as in Eq. (1.46). With 
growing N, the unbiased states increasingly differ from the states that we have 
used for the nonlocality test [see Eq. (1.47) for large squeezing]. On the other 
hand, the biased and the unbiased states are equivalent under local squeezing 
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operations and thus cannot differ in their potential nonlocality. In addition, 
this equivalence shows that also the unbiased states are only nonmaximally 
entangled for finite squeezing, which suggests that they also do not lead to an 
exponential increase of the violations as for the qubit GHZ states. 

In section 2.1, we discussed some properties of pure, fully entangled states 
of three qubits. An important feature of these states is that a distinction can 
be made between two inequivalent classes: states from the first class can be 
converted into the state |GHZ) via SLOCC, but not into the state | W) (not even 
with arbitrarily small probability). For the second class, exactly the opposite 
holds. In several senses, the representative |GHZ) of the former class would be 
best described as a maximally entangled state, whereas the representative |W) 
of the latter class is nonmaximally entangled. A distinct feature of the maximum 
entanglement of |GHZ) is that after tracing out one qubit, the remaining qubit 
pair is in a separable mixed state 8 . Apparently, the entanglement of |GHZ) 
heavily relies on all three parties. By contrast, the entanglement of the state 
|W) is robust to some extent against disposal of one qubit. When tracing 
out one qubit of |W), the remaining pair shares a mixed entangled state. In 
the continuous-variable setting, we can make analogous observations. By 
interpreting the state f dx\x, x, x) as the analogue of |GHZ), we see that 



Tri / dx dx' \x, x, x)(x', x' , x'\ = \ dx \x)22(x\ <8> |s)33(a;| , (1-79) 



which is clearly a separable mixed state (and indeed not the maximally mixed 
state oc J dx dx' \x, x') (x, x'\). More interesting is the behaviour of a regular- 
ized version of f dx \x, x, x). In order to apply bipartite inseparability criteria, 
let us trace out (integrate out) one mode of the Wigner function H / out (x, p) in 
Eq. (1.50) for N = 3, 




) 




oc cxp - 2e +2r 



e +2r + 2e -2r 

e" 2r + 2e+ 2r 



(4 + s§) - 2e~ 2r 



e -2r + 2e +2r 

e+ 2r + 2e~ 2r 




e 



+2r 



e 



e ~2r _ e +2r 



+4e 



,+2r 



x 2 x 3 + 4e 



P2P3 



e 



+ 2e+ 2r 



e+ 2r + 2e" 2r 



(1.80) 



"However, remember that this is not the maximally mixed state for two qubits. Only when tracing out two 
parties do we end up having the maximally mixed one-qubit state. 
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From the resulting Gaussian two-mode Wigner function, we can extract the 
inverse correlation matrix. For Gaussian iV-mode states with zero mean values, 
the Wigner function is given by 

H/(e= (^4a5v exp {4^ T }- <U1 > 

with the 2iV-dimensional vector £ having the quadrature pairs of all N modes 
as its components, 

£ = (x 1 ,p 1 ,X 2 ,P2, ...,X N ,p N ) , £ = (xi,P!,X2,P2, ...,X N ,p N ) , (1.82) 

and with the 2N x 2N correlation matrix V having as its elements the second 
moments (symmetrized according to the Weyl correspondence), 

TY[p(A£a4- + a4-A£)/2] = {(iiij + ijii)^) 

W(£)^ 3 d 2N t = V ij , (1.83) 



/ 



where A£j = £j — = & for zero mean values. The last equality defines 
the correlation matrix for any quantum state, but for Gaussian states of the 
form Eq. (1.81), the Wigner function is completely determined by the second- 
moment correlation matrix. Now we can calculate the bipartite correlation 
matrix of the state in Eq. (1.80), 

/ e +2r + 2e~ 2r 2 sinh 2r \ 

e- 2r + 2e +2r -2sinh2r 

2sinh2r e +2r + 2e~ 2r 

^ -2sinh2r e~ 2r + 2e +2r J 

(1.84) 

We could have also obtained this two-mode correlation matrix by extracting 
the three-mode correlation matrix V of the state VFo Ut (x, p) in Eq. (1.50) with 
N = 3 and ignoring all entries involving mode 1 [or equivalently by explicitly 
calculating the correlations between modes 2 and 3 with the Heisenberg oper- 
ators in Eq. (1.18) for r = r\ = r 2 and iV = 3]. The resulting two-mode state 
is a (mixed) inseparable state for any nonzero squeezing r > 0. Note that, for 
instance, the total variance in Eq. (1.37) with N = 2 becomes for this state 
(5e~ 2r + e +2r )/6, which drops below the boundary of 1 only for sufficiently 
small nonzero squeezing, but approaches infinity as the squeezing increases. 
However, we can easily verify the state's inseparability for any r > by look- 
ing at the necessary two-party separability condition in product form given in 
Ref. [37]. We find that 

([A(x 2 - x 3 )] 2 )([A(p 2 +p 3 )} 2 ) = (2e- 4r + 1)/12 , (1.85) 
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which drops below the separability boundary of 1/4 for any r > 0. Of course, 
also the necessary and sufficient partial transpose criterion from Ref. [14] 
indicates entanglement for any r > [28]. Recall that by first taking the 
"infinite-squeezing limit" and then tracing out one mode, we had obtained a 
separable state [Eq. (1.79)]. That was what we expected according to the result 
for the maximally entangled qubit state |GHZ). 

So after all, we confirm what we had intuitively expected: the tripartite 
state Wout( x ) p) far finite squeezing is a nonmaximally entangled state like the 
qubit state |W). Only for infinite squeezing does it approach the maximally 
entangled state J dx\x,x,x), the analogue of |GHZ). This result reflects 
what is known for two parties. The two-mode squeezed state W ou t( x 5 p) with 
N = 2 becomes a maximally entangled state J dx\x,x), such as the Bell state 
(|00) + | ll))/\/2, only for infinite squeezing. For finite squeezing, it is known 
to be nonmaximally entangled. 

3. CONCLUSIONS 

Let us conclude by asking whether we were able to find answers to the 
questions posed at the beginning of this chapter: how can we generate, measure, 
and (theoretically and experimentally) verify genuine multipartite entangled 
states for continuous variables? How do the continuous- variable states compare 
to their qubit counterparts with respect to various properties? 

As for the generation, we demonstrated that genuinely iV-party entangled 
states are producible with squeezed light resources and beam splitters. In 
particular, one sufficiently squeezed light mode is in principle the only resource 
needed to create any degree of genuine multi-party entanglement by means of 
linear optics. The resulting states, though genuinely multi-party entangled, are 
always nonmaximally entangled multi-party states due to the finite amount of 
the squeezing. They behave like the iV-party versions of the qubit state |W). 
First, they also contain bipartite entanglement readily available between any 
pair of modes (just as |W) and as opposed to the qubit state |GHZ)). Secondly, 
they yield a non-exponential increase of violations of multi-party Bell-type 
inequalities for growing number of parties (as for |W) and different from the 
qubit state |GHZ) for which the increase is exponential). 

Furthermore, we have seen that by inverting the circuits for generating 
genuine but nonmaximum multi-party entanglement, one can perform projec- 
tion measurements onto the maximally entangled multi-party (GHZ) basis for 
continuous variables. In contrast to the difficulties in performing such measure- 
ments for photonic qubits within the framework of linear optics, continuous- 
variable GHZ measurements only require beam splitters and homodyne detec- 
tors. In addition, we showed that the circuits for measuring maximum GHZ 
entanglement are also applicable to the theoretical and experimental verifica- 
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tion of the nonmaximum entanglement of the multi-party states (which are 
those producible in the laboratory). The circuits provide a necessary condition 
for full separability of any iV-partite iV-mode state (pure or mixed, Gaussian 
or non-Gaussian) with any number of modes N. However, this condition is 
not sufficient for full separability and, more importantly, its violation does not 
verify genuine but only partial multipartite entanglement. For the theoretical 
verification of genuine multipartite entanglement, additional assumptions have 
to be taken into account such as the total symmetry of the relevant states. There- 
fore, an unambiguous experimental proof of genuine multipartite entanglement 
of continuous-variable states was not proposed in this chapter. A possible ap- 
proach to this would be to consider the violation of stricter iV-party Bell-type 
inequalities which cannot be violated by only partially entangled states. How- 
ever, the experimental nonlocality test would then rely on observables such as 
the photon number parity, and hence become unfeasible with current technol- 
ogy. More desirable would be a test for genuine multipartite entanglement that 
is solely based on linear optics and efficient homodyne detections. 
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